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ABSTRACT 

An exceptional 7 -ray outburst from 3C 279 is detected by Lermi-Large Area Telescope (LAT) in 2015 June. 

In the energy range of 0.1—300 GeV, the highest flux measured is (39.1±2.5) x 10“® ph cm“^ s“^, which is 
the highest 7 -ray flux ever detected from 3C 279, exceeding the previous historically brightest flare observed 
by EGRET in 1996. The high activity period consists of three major flares with the last one being the brightest. 

All but one flares show a faster rise and slower decay pattern and at the peak of the activity, the 7 -ray spectrum 
is found to show a clear signature of break/curvature. The obtained spectral parameters hint for the peak of the 
inverse Compton emission to lie in the LAT energy range (around ~1 GeV) which is in contrast to that seen 
during the 2013 December and 2014 April 7 -ray flares of 3C 279. From the 77 pair opacity arguments, the 
minimum Doppler factor is estimated to be 14 and the location of the 7 -ray emitting region is found to be either 
at the outer edge of the broad line region or farther out from it. 

Subject headings: galaxies; active — gamma rays: galaxies — quasars: individual (3C 279) — galaxies: jets 


1. INTRODUCTION 

The radio source 3C 279 (z=0.536; Lynds et al. 1965) is 
one of the best studied flat spectrum radio quasars (FSRQ) 
in all the wavebands. It emits strong and variable emission 
at all frequencies and is one of the first blazars detected by 
EGRET instrument on the Compton Gamma-Ray Observatory 
(Hartman et al. 1992). Since 2008, 3C 279 is being continu¬ 
ously monitored by Eermi-Large Area Telescope (LAT) and 
was a subject of various multi-wavelength campaigns (e.g. 
Hayashidaet al. 2012, 2015). This source is an active 7 -ray 
emitter and multiple episodes of 7 -ray outbursts have been 
detected (Hayashida et al. 2012; Paliya et al. 2015). In partic¬ 
ular, two prominent flares of similar intensity (LceV ^ 12 ; 
where LbeV is in 0.1—300 GeV range and in units of 10“® 
ph cm“^ s“^) have been detected by LAT in 2013 Decem¬ 
ber and 2014 April. Interestingly, the 2013 December flare 
exhibited a hard rising spectrum, whereas a significant curva¬ 
ture was noticed in 2014 April flare (Hayashida et al. 2015; 
Paliya et al. 2015). A 7 -ray flux variability as short as ~1 hr 
is also reported (Paliya et al. 2015). 

In 2015 June, 3C 279 displayed an exceptionally high activ¬ 
ity state. The daily binned 7 -ray flux is reported as the highest 
measured since the beginning of the Fermi operation (Cutini 
2015). Two target of opportunity observations were approved 
for a total of 750 ksec duration (between MJD 57,188 and 
57,197) and during this period Fermi observed 3C 279 in 
a special pointed mode, other than its normal sky scanning 
mode operation. In this letter, motivated by the availability of 
good quality LAT data, the intraday 7 -ray variability and as¬ 
sociated spectral changes are studied in detail. The obtained 
results are further compared with the previous major flares de¬ 
tected from 3C 279. In Section 2, observations and data anal¬ 
ysis of Fermi-LKL data from 2015 June 7 to June 20 (MJD 
57,180—57,193) are reported. Results are presented in Sec¬ 
tion 3 and discussed in Section 4. Throughout, a ACDM cos¬ 
mology with the Hubble constant Hq = 71 km s“^ Mpc“^, 
Lljn = 0.27, and Ua = 0.73 is adopted. 
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2. re7?M/-LARGE AREA TELESCOPE OBSERVATIONS 

To analyze the LAT data (P7REP) covering the period of 7 - 
ray outburst, the standard data reduction procedure is adopted 
and here it is described in brief. The unbinned likelihood 
method included in the pylikelihood library of Science 
Tools (v9r33p0) is used to extract the information from 
the events belonging to the SOURCE class and energy range 
0.1—300 GeV and lying within 10° region of interest (ROI) 
centered at the position of 3C 279. A maximum likelihood 
(ML) test statistic TS=2Alog(£), where C represents the like¬ 
lihood function between models with and without a point 
source at the position of source of interest, is calculated to 
determine the significance of the 7 -ray signal. All the sources 
lying within the ROI and defined in the third Fermi-LPIL cat¬ 
alog (3FGL; Acero et al. 2015) are included and their spec¬ 
tral parameters are left free to vary during the model fitting. 
Sources lying within 10° to 15° from the center of the ROI 
are also appeared in the model file and their parameters are 
kept fixed to the 3FGL catalog values. A first run of the ML 
analysis is performed over the period of interest and all the 
sources with TS<25 are removed from further analysis. 

To characterize the variability properties of the source, light 
curves are generated using various time binnings (3 hr, 6 hr, 
12 hr, and 1 day) and over different energy intervals (0.1—300 
GeV, 0.1 — 1 GeV, and 1—300 GeV). Though the 7 -ray spec¬ 
tral shape of 3C 279 exhibit a definite curvature, all the light 
curves are produced by adopting a simple power law (PL) in 
each time bin, since the statistical uncertainties on the PL in¬ 
dices are smaller compared to those obtained from complex 
model fits such as broken power law (BPL). 

Various models are applied to fit the 7 -ray spectra of 3C 
279 and that includes a BPL (N{E) = Ao(i7/i7break)~^S 

with i = 1 if E < inbreak and i = 2 if E > inbreak), a 
log-parabola (N{E) = No{E/Ep)~°'~^ ^os{E/Ep)^ where a 
is the photon index at Ep, /3 is the curvature index and Ep 
is fixed at 300 MeV), and a PL model over logarithmically 
equally spaced energy bins with L kept fixed to the value fitted 
over the whole energy range. 


2 


Vaidehi S. Paliya 


All the time/energy bins with AF,y/Fj >0.5, where AF^ 
is the error estimate in the flux F^, and/or TS<9' are rejected 
from the analysis. Systematics on the measured fluxes are of 
around 10% below 100 MeV, 5% between 316 MeV to 10 
GeV and 10% above 10 GeV. The statistical uncertainties are 
estimated at Icr level. 


the fitting is performed by including a constant background 
and three fast temporally evolving components. Each flare is 
approximated to be fit by a function of the following form 


F{t) = 2Fo 


exp 
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3. RESULTS 

The daily binned light curve of 3C 279 at photon energies 
0.1 —300 GeV, covering the period of high activity is shown in 
the top panel of Figure 1 . The source started showing activity 
around MJD 57,184 with flux level about FceV ~ 1 and after a 
period of 3 days, 7 -ray flux increased by a factor of 10 where 
it remain steady for 2 days. The source displayed another flare 
of higher amplitude on MJD 57,189 when the highest 7 -ray 
flux obtained as FceV = 24.5±0.5, which is the highest daily 
binned 7 -ray flux measurement from 3C 279 since the launch 
of Fermi-LAT. The associated photon index is hard and having 
a value of 2.05±0.02. Immediately after this flaring activity, 
the source returned to a relatively low activity state where flux 
remained at the level of FceV '^T Based on the observed 
activity in the daily binned light curve, the entire period is 
then divided into four states: a pre—flare, two flares, and a 
post—flare state. These periods are shown in the top panel of 
Figure 1 and labeled as Pre, PI, P2, and Post, respectively. 

Due to good photon statistics, the light curves are also gen¬ 
erated using time bins of 12 hr, 6 hr, and 3 hr and are shown in 
the bottom three panels of Figure 1. As can be seen, the flare 
of period PI starts getting resolved in finer bin light curves. 
In fact, this period comprises of two flares distinctly visible 
in 6 hr binned light curve. There is a third sharp flare as PI 
ends and P2 begins. These three flares are named as FI, F2, 
and F3 for the rest of the analysis (see Figure 1). The rising 
phase of the brightest flare F3 is clearly unresolved down to 
3 hr scale whereas its decaying segment seems to be nicely 
resolved in 6 hr binning. Though the flares FI and F2 ap¬ 
pear to be resolved in 6 hr binning, the 3 hr binning light 
curve suggests for the presence of sub-structures consisting of 
several large amplitude yet unresolved and exhibiting chaotic 
events. The highest 7 -ray flux, using 3 hr binning, is ob¬ 
tained as FoeV = 33.5±1.4. Moreover, the data is also an¬ 
alyzed using the time bins defined as Good Time Intervals 
(GTI, Foschini et al. 201 1). The highest 7 -ray flux measured 
using this method is FceV = 39.1 ±2.5, thereby making this 
peculiar 7 -ray flare as the brightest event recorded from 3C 
279 not only since the beginning of Fermi operation but also 
ever. For a comparison, the historically brightest flare was 
detected by EGRET when the measured flux was as high as 
Ki2 X 10“® ph cm“^ s“^(Wehrle et al. 1998). Further, the 
associated GTI bin size is ~14 min with TS = 2407 and 275 
counts are registered. 

To determine the duration of the shortest flux variability, the 
light curves are scanned using the following equation 

F{t2) = ( 1 ) 

where F[ti) and F{t 2 ) are the fluxes at time ti and t 2 re¬ 
spectively, and Td is the flux doubling/halving timescale. The 
shortest flux doubling time using this method is measured as 
2.2±0.3 hr on MJD 57,189 with ^9cr significance. Moreover, 
in order to assess the asymmetry of rise and decay of flares, 
all the three flares (FI, F2, and F3) are also subjected to tem¬ 
poral profile fitting. This is done for 6 hr binned flares and 

* TS of 9 corresponds to detection (Mattox et al. 1996). 


(Abdo et al. 2010), where Fq is the flux at time to representing 
approximately the flare amplitude, and and Td are the rise 
and decay time of the flare. The obtained rising time for FI, 
F2, and F3 are 5.40F0.84 hr, 4.05F2.73 hrs, and 3.06F0.57 
hr respectively, whereas the associated decay time of the flares 
are 9.91F5.04 hr, 4.13F2.02 hr, and 8.68F0.42 hr respec¬ 
tively. The result of the time profile fitting is shown in Fig¬ 
ure 2 . 

FAT 7 -ray light curves of 3C 279, focusing on the flare pe¬ 
riod, in two different energy bands: 0.1 — 1 GeV and 1—300 
GeV are shown in the top panel of Figure 3. The 1—300 
GeV light curve is scaled appropriately to compare its vari¬ 
ability pattern with that seen in 0.1 —1 GeV range. As can be 
seen, though the pattern of flux variations are similar in both 
the bands during F3 flare, there are very moderate change in 
source brightness in 1—300 GeV band compared to 0.1 — 1 
GeV for FI and F2 events. This observation indicates that 
the flare FI and F2 are primarily caused by low energy elec¬ 
trons. Further, the hardness ratio is also calculated to assess 
the spectral evolution as a function of time. This is done using 
following equation 


.Fh — Fs 




Fh ± Fs 


(3) 


where Fs and Fh are 6 hr binned 7 -ray fluxes in 0.1 — 1 GeV 
and 1—300 GeV energy ranges, respectively. Though there 
are no significant spectral hardening/softening during FI and 
F2 flares, the spectrum is clearly hard at the peak of F3 and 
it softens as the flare decays (see middle panel of Figure 3). 
Further, recently it has been proposed that the origin of 7 -ray 
flares to be lying inside/outside the BLR is reflected from the 
absence/presence of time lag between MeV and GeV emis¬ 
sion (Dotson et al. 2012). To test this hypothesis, a time lag 
analysis is performed using the z-transformed discrete corre¬ 
lation function (ZDCF) method of Alexander (2013) (see also 
Alexander 1997). The errors are computed using a Monte 
Carlo simulation by adding a random error at each step to each 
data from the errors in the light curves (Alexander 2013). The 
result of this analysis is shown in bottom panel of Figure 3 . As 
can be seen, there is no significant lead/lag noticed between 
two energy bands (time lag of O.O^q days). It should be 
noted that there could be lag between these light curves, how¬ 
ever, limited time resolution of the light curves and closeness 
of the FAT energy bands makes it difficult to quantify. 

The hints for the presence of spectral evolution can be fur¬ 
ther investigated by plotting photon index versus flux. This 
is done for 6 hr binned data, in two energy ranges: 0.1—300 
GeV and 0.1 — 1 GeV, for the period PI and P2 (covering the 
flares FI, F2, and F3). The results are presented in the left 
panel of Figure 4. For comparison, the averages obtained 
over the period ‘Pre’, is also shown with blue squares. A 
clockwise pattern is evident during P2 period for both energy 
bands. Though there are hints for the same behavior during 
PI period in 0.1 —300 GeV energy range, the 0.1 — 1 GeV plot 
show a clear evidence of clockwise helical pattern. The re¬ 
duced for a constant fit of the photon index for the period 
PI are 6.94/7 and 6.68/7, for 0.1-300 GeV and 0.1-1 GeV 
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respectively. Comparing to ‘Pre’ period, a flux enhancement 
by a factor of 10 is accompanied by slight hardening of the 
spectrum. While drawing any firm conclusion on particle ac¬ 
celeration and cooling from these patterns is not possible, the 
lack of strong spectral variability still provides clues to the 
underlying radiative processes. Due to poor statistics, a clear 
pattern of the photon index for energy range 1—300 GeV ver¬ 
sus flux could not be observed with 6 hr time binning. 

In the right panel of Figure 4, the LAT spectra obtained 
by averaging the four periods shown in Figure 1, are dis¬ 
played. These flux distributions have been fitted with PL (dot¬ 
ted), log-parabola (dashed), and BPL (solid) models. Due 
to low photon statistics during the periods ‘Pre’ and ‘Post’, 
the BPL model could not be fitted. To substantiate the pres¬ 
ence/absence of the spectral curvature, the TS of the curvature 
TScnrve = 2(log £(log-parabola/BPL) — log £(PL)) is also 
computed (Nolan et al. 2012). The associated parameters are 
provided in Table 1 . Out of all the four periods considered, 
a statistical significant break/curvature is noticed only during 
PI and P2 periods. Both the log-parabola and BPL model re¬ 
produce the data satisfactorily, though the log-parabola model 
gives slightly better fit. The break energy obtained from the 
BPL model fitting remains constant irrespective of the fluctu¬ 
ations seen in the flux levels, during both the periods. This 
finding is inline with that observed during the 7 -ray outburst 
of 3C 454.3 in 2010 November (Abdo et al. 201 1). 


4. DISCUSSION 


The 7 -ray outburst of 3C 279 in 2015 June has broken 
the record of the highest 7 -ray flux ever measured from this 
blazar. The entire flaring episode is found to comprise of three 
flares with the last one having the largest amplitude. The ris¬ 
ing part of the brightest flare is unresolved down to 3 hr bin¬ 
ning and this hints that with the limited time resolution, the 
apparent profiles of the source flux variations may not reflect 
the exact temporal characteristics of the 7 -ray outburst. How¬ 
ever, the decaying part of the flare is nicely resolved in 6 hr 
and 3 hr binned light curves. 

The shortest flux doubling timescale obtained in this work 
is ~2 hr which is similar to that obtained during 2013 De¬ 
cember outburst of 3C 279 (Hayashida et al. 2015). Interest¬ 
ingly, comparing to 2014 April outburst, though the flux level 
is significantly higher during 2015 June flare, the faster vari¬ 
ability was measured during the former (Paliyaet al. 2015). 
However, similar to 2013 December and 2014 April flares, the 
flares of 2015 June outburst are found to exhibit asymmetric 
variability patterns with fast rise and slow decay (except for 
F2), which can be explained by invoking the rapid injection of 
accelerated electrons, probably at shock front, and the decay 
can be attributed to the weakening of the shock. 

The minimum Doppler factor (jmin can be estimated numer¬ 
ically from 77 opacity arguments and by measuring the en¬ 
ergy of the highest energy photon detected during the flare. 
Assuming that the optical depth T^j{ei) of a photon with en¬ 
ergy ei = Ei/meC^ to the 77 interaction is = 1 , the min¬ 
imum Doppler factor can be calculated as follows (see e.g. 
Dondi & Ghisellini 1995; Ackermann et al. 2010) 




cTTd\{l -f 




rlTleC* 


1/6 


(4) 


where /e is the flux at energy e = hvjmeC^ and e = e = 
25^/(l + z)^ei (Ackermann et al. 2010). In this calculation, 
the energy of the highest energy photon (Ei) is used during 


the period when variability time fvar and are measured. 
The highest energy photon of energy Ri52 GeV is detected 
on MJD 57,189.62 at an angular separation of 0°.09 from the 
3FGL position of 3C 279, with 99.99% probability of being 
source photon (see bold downward arrow in the bottom panel 
of Figure 1). Almost at the same time, the Swift X-ray Tele¬ 
scope observation has revealed the energy flux in 0.3—10 keV 
band as 3.5 x 10~^^ erg cm“^ and the associated pho¬ 
ton index is 1.4 (Pittori et al. 2015). Since the highest energy 
photon was detected during the fading part of F3, the variabil¬ 
ity time can be taken as fvar = ln(2) « 6 hr. This gives 

« 14. 

A rough estimation about the location of the emission re¬ 
gion (i?diss) can be done by assuming it to have spherical 
geometry and bulk Lorentz factor F = Jmin = 14. Then 
i?diss < 2cr^fvar/(l + z:) Ri0.05 pc. Now, Considering 
the accretion disk luminosity (Ldisk) of 3C 279 as 1 x 10^® 
erg s“^ (Paliya et al. 2015), the size of the BLR comes out to 
be ~1 X 10“ cm Ri0.03 pc (e.g. Kaspi et al. 2007). There¬ 
fore, at the time of the 52 GeV photon emission, the 7 -ray 
radiating region must have been located either close to the 
outer edge of the BLR or beyond it. This is also supported 
by fact that to avoid severe absorption of high energy photons 
(having energy >50 GeV) by optical-UV radiation of BLR 
radiation field, the emitting region has to be farther out from 
the BLR. 

The highest 7 -ray flux measured is EceV = 39.1±2.5 and 
the associated photon index is 2.0±0.L This corresponds to 
an isotropic luminosity (Lj) of (5.5±L1) x 10"^® ergs“^. 
Assuming F = 14, the 7 -ray luminosity in the jet frame 
would be ~ 1.4 x 10'^^ ergs“^. This is 

about 10 times the total available accretion power (Lacc — 
J^disk/i/disk — 1 X 10^® erg s“^; assuming radiative efficiency 
? 7 disk = 10%). The black hole mass of 3C 279 is found in 
the range of (3—8) x 10 ®Mq (Gu et al. 2001; Woo & Urry 
2002), and thus the Eddington luminosity LEdd ~ (0.4—1) 
X 1047 ergs 4. Now, in order for the observed to be 
lower than LEdd, 0, < 5“ - 7° and a beaming factor (1 — 
cos 0j)“4 > 240 — 440. Similar parameters are reported by 
Abdoet al. (2011) for the GeV flare of 3C 454.3 during its 
2010 November flare. 

The cooling timescale of the 7 -ray (e-y = 1 GeV) emitting 
electrons can be calculated as follows (e.g. Saito et al. 2013) 


SnipC 


A;ool — 


' co(l + z) 


4cttU(: 


BLR/torus 


(5) 


where uj^LR/torus = /BLR/torusidiskrV 47 ri? 2 ^^/^^^^^, is 
the comoving seed photon energy density for the BLR/torus 
radiation field with characteristic energy eq = 10.2/0.2 eV and 
ZeLR/torus is the fraction of Ldisk reprocessed by BLR/torus. 
This resulted in TcooI = 7 min or 11 min, provided the 
seed photons for EC process are originated from the BLR 
or torus, respectively. The obtained cooling time is signifi¬ 
cantly shorter than that observed decay timescale of the flares, 
suggesting that the observed flare decrease is controlled not 
only by radiative cooling, but by a combination of various 
factors such as the geometry and sub-structure of the emis¬ 
sion region or possibly the jet dynamics (Tanihata et al. 2001; 
Kushwaha et al. 2014). 

A significant break/curvature in the 7 -ray spectrum is the 
characteristic properties of powerful ESRQs. Interestingly, 
this feature is found to be more prominent during the flaring 
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Table 1 

Results of the Model Fitting to the 7 -ray Spectra of 3C 279, obtained for different activity periods. 




Power law 





Activity 

FOeV 

To. 1-300 GeV 



TS 



( 10 “® ph cm“^ s“l) 






Pre 

1.41±0.19 

2.27±0.11 

- 

- 

274.55 

- 

PI 

15.80±0.55 

2.21±0.03 

- 

- 

6553.42 

- 

P2 

24.50±0.49 

2.05±0.02 

- 

- 

22673.92 

- 

Post 

2.39±0.13 

2.26±0.05 

- 

- 

1690.61 

- 



Log parabola 





Activity 

FoeV 

a 

P 


TS 

T (Scurve 


( 10 “® ph cm~^ s“l) 






Pre 

1.41±0.19 

2.27±0.11 

O.OOiO.OO 

- 

274.55 

0.00 

PI 

15.20±0.55 

2.05±0.05 

0.11 ±0.03 

- 

6572.25 

18.83 

P2 

23.60±0.48 

1.84±0.03 

0.13±0.02 

- 

22769.93 

96.01 

Post 

2.34±0.13 

2.16±0.08 

0.08F0.04 

- 

1694.40 

3.79 



Broken power law 





Activity 

FoeV 

Ti 

T2 

-^break 

TS 

T S'curve 


( 10 “® ph cm“^ s“l) 



(GeV) 



PI 

15.31±0.55 

2.05±0.05 

2.60F0.11 

l.UU_o 30 

6571.30 

17.88 

P2 

23.85±0.48 

1.88±0.03 

2.58F0.07 

+-4^0.i5 

22757.23 

83.31 


episodes (see e.g. Abdo et al. 2011; Paliya et al. 2015). Dur¬ 
ing P2 period, i.e. at the peak of the flare, the BPL photon 
index Fi before inbreak (~1 GeV) is < 2, thus indicating a 
rising spectrum, whereas ^2 > 2, implying a falling spectrum 
after inbreak- This observation suggests that probably the peak 
of the IC mechanism is seen in the LAT bandpass and thus 
the observed shape most likely reflects the energy distribu¬ 
tion of the emitting electrons. Comparing the shape of the 7 - 
ray spectrum with that seen during 2013 December and 2014 
April reveals an interesting observation. The 2013 Decem¬ 
ber flare exhibited a hard rising spectrum (F = 1.71 ± 0.10, 
Hayashida et al. 2015) indicating that the IC mechanism to 
peak at very high frequencies, whereas 2014 April events 
showed a steep falling spectrum (F = 2.23±0.03, Paliya et al. 
2015), thus implying IC peak to lie before the LAT energy 
range. These three flares, therefore, represent the three dif¬ 
ferent characteristics of the radiative processes powering the 
relativistic jet of 3C 279. The observation from Fenni-LAJ', 
in this regard, is an invaluable asset to learn the physical prop¬ 
erties of the most powerful objects in the Universe, called 
blazars. 
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Figure 2. Temporal profile fitting of 6 hr binned 7 -ray light curve. Thi'ee flares are shown by cyan, pink, and limegreen colors, whereas grey line represnts the 
constant background. Black line is the sum of all the components. 
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Figure 3. Top: Six hr binned light curves of the fluxes in 0.1 — 1 GeV (red) and 1—300 GeV, multiplied by a factor of 12 (black). Middle panel presents the 
variation of HR. Bottom: DCF calculated for 0.1 — 1 GeV and 1—300 GeV light curves. The time ordering is T(( 1—300 GeV)—(0.1 —1 GeV)). 
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Figure 4. Left: T vs. flux for the periods PI and P2, obtained with 6 hr binning for 0.1—300 GeV and 0.1 — 1 GeV bands. Right: Fermi-LA^ SEDs during 
different activity states as defined in Figure 1. PL, LP, and BPL models are shown with dotted, dashed, and solid lines respectively. Horizontal errorbar 
corresponds to energy range of each bin whereas vertical bar represents Icr statistical errors. The residuals in the lower panel refer to the PL model. 













































